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Abstract 

We study random domino tilings of a Double Aztec diamond, a 
region consisting of two overlapping Aztec diamonds. The random 
tilings give rise to two discrete determinantal point processes called the 
K- and L-particle processes. The correlation kernel of the K-particles 
was derived in Adler, Johansson and van Moerbeke (2011), who used 
it to study the limit process of the K-particles with different weights 
for horizontal and vertical dominos. Let the size of both, the Double 
Aztec diamond and the overlap, tend to infinity such that the two 
arctic ellipses just touch; then they show that the fluctuations of the 
K-particles near the tangency point tend to the tacnode process. In 
this paper, we find the limiting point process of the L-particles in the 
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overlap when the weights of the horizontal and vertical dominos are 
equal, or asymptotically equal, as the Double Aztec diamond grows, 
while keeping the overlap finite. In this case the two limiting arctic 
circles are tangent in the overlap and the behavior of the L-particles 
in the vicinity of the point of tangency can then be viewed as two 
colliding GUE-minor process, which we call the tacnode GUE minor 
process. As part of the derivation of the kernel for the L-particles 
we find the inverse Kasteleyn matrix for the dimer model version of 
Double Aztec diamond. 
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1 Introduction and main results 

The problem of a random domino tiling of a single Aztec diamond has been 
widely investigated by the combinatorics and probability community; the 
highlight was the existence of an inscribed arctic circle: inside the circle the 
dominos display a disordered pattern and outside a regular brick wall pattern; 
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2( 2n -jo + l) 



Figure 1: Double Aztec diamond with n = 8 and overlap p = 4, with the 
(£, 77) coordinates with Aztec diamond A enclosed by the blue dotted line 
and Aztec diamond B enclosed by the red dotted line. The overlap contains 
p lines (through black squares) £ = 2s for n — p < s < n. 
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see [H El EH El [151 [13 • When the weight a of vertical dominos is different 
from the one of horizontal dominos, then the arctic circle gets replaced by 
an inscribed arctic ellipse. 

In [2] the authors investigate the domino tiling of two overlapping Aztec 
diamonds, each of size n, with weight < a < 1 for vertical dominos and 
weight 1 for horizontal dominos. When the size of the diamonds and the 
overlap both become very large, in such a way that the two arctic ellipses 
for the single Aztec diamonds merely touch, then a new critical process, the 
tacnode process, will appear near the point of osculation (tacnode); it is run 
with a time in the direction of the common tangent to the ellipses. The kernel 
governing the local statistics of the tacnode process is given by a perturbation 
of the Airy process kernel by an integral of two functions. It was also shown 
in [2] that this tacnode process has some universal character: it coincides 
with the one found in the context of two groups of non-intersecting random 
walks [1] and Brownian motions, meeting momentarily [HI [10]; see also [5]. 

Another ingredient here is the process given by the successive interlacing 
eigenvalues of minors of a GUE-matrix: the so-called GUE-minor process. 
In [IT] this process has arisen in the following model: magnifying the in- 
finitesimal region about the point of tangency of the arctic ellipses with the 
edge of a single Aztec diamond for large n, leads to a determinantal process 
of interlacing points on the successive lines through (say) the black squares, 
parallel to the edge of the diamond. This yields the GUE-minor process, see 
also [22]. 

In the present work, we consider two overlapping Aztec diamonds with 
an overlap, which remains finite, when the size of the diamonds tends to oo. 
In order to maintain the osculation of the two inscribed ellipses, the geom- 
etry forces the weight a of the vertical ones to tend to the weight 1 of the 
horizontal (3/ y/n/2. Macroscopically this amounts to two 

Aztec diamonds with inscribed arctic circles intersecting infinitesimally. In 
view of the comments above, it seems natural that this process be related 
to the GUE-minor kernel. Indeed, when n — > oo, looking with a magni- 
fying glass at the infinitesimal overlap of the diamonds gives rise to a new 
determinantal process, with local statistics given by the so-called tacnode 
GUE-minor kernel; it is a finite rank perturbation of the GUE-minor kernel 
mentioned above. As far as we are aware, there is no Random Matrix theory 
counterpart of this distribution. 

In [2], the authors considered successive lines through black squares per- 
pendicular to the region of overlap with dots in the black square each time 
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the dominos covering that square is pointing to the right of or above the 
line; these are called the K-particles. In this work, we shall mainly consider 
successive lines parallel to the region of overlap and put a dot in the black 
square each time the dominos covering that square is pointing to the left of 
or above that line; these give the L-particles (Section 1.2). We deduce the 
kernel for the L-particles from the one of the K-particles, by first obtaining 
the inverse Kasteleyn matrix [19] for the double Aztec diamond in terms of 
the K-kernel and then deduce the L-kernel of the particles from that same 
inverse Kasteleyn matrix (Section 2). The inverse Kasteleyn matrix of a sin- 
gle Aztec diamond had been obtained for a = 1 by [IB] and generalized for 
all a recently in [3]. 

In Section 1.2, we study the specific interlacing pattern of the L-particles. 
We state in Section 1.3 and prove in Section 5 that, in the scaling limit 
n — > oo, the L-process in the infinitesimal overlap is indeed driven by the 
tacnode GUE-minor kernel. Also, we state in Section 1.3 and prove in Section 
5 that upon thinning at the rate p n = 1—2/ the K-process is also driven 
by the same tacnode GUE-minor kernel, but with a (somewhat surprising) 
shift in one of the discrete parameters. 

1.1 Domino tilings of double Aztec diamonds and ran- 
dom surface 

Consider two overlapping Aztec diamonds A and B, of equal sizes n and 
overlap p, with opposite orientations; i.e., the upper-left square for diamond 
A is black and is white for diamond B. The size n is the number of squares 
on the upper-left side and the amount of overlap p counts the number of 
lines of black squares common to both diamonds A and B. Let £, rj be a 
system of coordinates as indicated in Figure [T] The even lines £ = 2k for 
< k < 2n — p and the odd lines T] = 2k — 1 for 1 < k < n run through 
black squares. The p even lines £ = 2{n — p + 1), ■ ■ ■ , 2n belong to the overlap 
of the two diamonds. Cover this double diamond randomly by dominos, 
horizontal and vertical ones, as in Figure [3j The position of a domino on 
the Aztec diamond corresponds to four different patterns, given in Figure [3] 
below: North, South, East, West. 



Define 




white squares of diamond A along the line rj = 
having an edge in common with the boundary 
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Figure 2: Double Aztec diamond with n = 8 and overlap p = 4 with the 
(z, x) coordinates of the Aztec diamond. 



and define m such that 



M- 1 



n — p 



2m 



2m 



when n and p have same parity 
— 1, when n and p have opposite parity. 



Throughout the paper we assume, for simplicity, the same parity for n and 
p, so that n — p = 2m. 

Together with this arbitrary domino-tiling of the double Aztec diamond 
A U B, one defines a piecewise-linear random surface, by means of a height 
function h specified by the heights, prescribed on the single dominos ac- 
cording to figure [4] above; this height can be taken to be piecewise-linear on 
each domino. This height function is different from the usual one by Cohn, 
Kenyon, Propp |4J, but related to it by an affine transformation. Let the 
upper- most edge of the double diamond A U B have height h = 0. Then, 
regardless of the covering by dominos, the height function along the bound- 
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North East South West 



Figure 3: Random Tiling of a double Aztec diamond with n = 8 and p — 4. 
The figure on the left shows the underlying checkerboard structure while the 
figure on the right shows the same tiling with the level lines. These are shown 
explicitly for the four types of dominos in the bottom figure. These level lines 
are the same as the DR lattice paths [21 J . 



ary of the double diamond will always be as indicated in Figure |4[ with 
height h = 2n along the lower-most edge of the double diamond. Away from 
the boundary the height function will, of course, depend on the tiling; the 
associated heights are given in Figure |4[ 

The height function h obtained in this way defines the domino tiling in 
a unique way, because a white square together with its height specifies in a 
unique way to which domino it belongs to: North, South, East and West; 
the same holds for black squares. 

This height function associates thus a piece-wise linear random surface 
with each random tiling and two groups of level curves of this random surface 
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2 2 
3 3 




h h h h + 1 (t+l/t+lA+lA+1 4+1 /; + 1 

North East South West 



Figure 4: The level lines including the height function around the boundary. 
The heights change in the interior only when crossing a level line. The bottom 
figure shows the height change for each individual domino. 
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corresponding to the half-integer values: 

13 5 11 1 

-, -, n , n-\ — 2n . 

2 2 2 2 2 2 

S v ' S v ' 

v4-level curves £>-level curves 

Put the weight a > on vertical dominoes and the weight 1 on horizontal 
dominoes, so that the probability of a tiling configuration T can be expressed 
as 

fl #vertical domino's in T 

Pfdomino tiling T) = -. — t—, : — -. — : 

fl #vertical domino s m t 

all possible tilings T 

Remember 2m := n — p throughout the paper. We will also use the coor- 
dinates indicated in figure [2} These are the coordinates, which we will call 
diamond coordinates, that were used for the particle processes in [2]. The 
transformation from diamond coordinates (z, x) to Kasteleyn coordinates 

(6, v) is s iven by : 

z = r] + 1 !; = z — 2x + 2m 

x = \{r] — £ + 2m + 1) rj = z — 1 



1.2 Two determinant al point processes L and IK 

I. The h-process is specified by putting a dot in the middle of the black 
square when the line £ = 2s in (£, -^-coordinates for < s < 2n — p intersects 
a level curve. We call these dots h-particles. See Figure [5] for an example. 
More precisely we can put a blue dot when intersecting A-level curves and 
a red dot when intersecting 5-level curves to distinguish the dots coming 
from the two Aztec diamonds; see Figure [6] In other terms, put a dot in the 
black square each time the random surface goes down one unit along the line 
£ = 2s. We are concerned with the probabilities of the following kinds of 
events, where [k, £] is an interval of odd integers along the maxis (so k and i 
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can be taken odd): 



{The line {£ = 2s} has an i]-gap D [k,£]} 

= {Interval [k,£] C {£ = 2s} in ^-coordinates contains no dot-particles} 
= {The random surface is flat along the ^-interval [k,£] C {£ = 2s}} 
= {Dominos covering [k,£] C {£ = 2s} are pointing to the left of 
or above the line {£ = 2s}} 
= {Dominos covering [k,£] C {£ = 2s} are red or yellow in upper Figure [5 



Theorem 1.1 The h-particles on the successive lines {£ = 2s} for 1 < s < 
2n — p form a determinants point process with correlation kernel 

L n , p (a, w 6, m) =(i + « 2 )4 0) (6, tk 6, m) 

- (1 + a 2 )(((J - K n )-l p+i A^ m )(k), B^ m (k))^ p+1 . 

given by a perturbation of a kernel by an inner-produc^ involving the 
ilvent of yet another kernel K n , all given by formulas (section {2.1 )) 



resoL 



This shows that, given q lines {£ = 2sj} and integers k i} £ i; with < k{ < 
li < n, the gap probability is expressed as the Fredholm determinant 

P | C\ {The line {£ = 2sj} has an r\-gap D 

\i=l / 

= det (l - [xM](^)W( 2s <.^; 2s j>%-)X[fc i A]('?i)]i<y<, 
of the kernel L n)P . 

This theorem will be proved in section |2j 

II. The WL-process. Now we put instead a blue dot in the middle of 
the black square when the line z = 2k in (z, x)-coordinates for 1 < k < n 
intersects an A-level curve and a red dot when intersecting a S-level curve; 
i.e., put a dot each time the random surface goes down one unit along the line 
z = 2k; see Figure [5] These dots define the K-particles. In this instance, we 
are concerned with the probabilities of the following kinds of events, where 

1 (/( fc ))5(fc))>a = f( k )g( k ) is an inner product in £ 2 {a, oo]. 
2 The variables r\i below run through odd values only. 
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Figure 5: The top figures show the L particle process using the level lines (on 
the left) and the dominos on the right. There is an L particle for every green 
(south) and blue (east) domino. The bottom figures show the IK particle 
process. There is a K particle for every green (south) and red (west) domino. 
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[k, £} is an interval along the x-axis: 
{The line {z = 2r} has an x-gap D [k, £]} 

= {The interval [k,£] C {z = 2r} in x-coordinates contains no dot-particles} 
= {The random surface is flat along the x-interval [k,£] C {z = 2r}} 
= {Dominos covering [k,£] C {z = 2r} are pointing to the left 
or below {z = 2r}} 

= {Dominos covering [k,£] C {z = 2r} are blue or yellow in lower Figure [5]f 



Theorem 1.2 [2 J The WL-particles on the successive lines {z = 2r} for 1 < 
r < n form a determinantal point process with correlation kernel given by 
perturbing the one-Aztec diamond kernel K° with an inner-product involving 
the resolvent of the kernel K n , all defined in {2.1 ) and |5p: 

(-l) x - y K n , p (2r,x;2s,y) = K° (2r, x; 2s, y) 

> ' >n-p+l 

This shows that given q lines {z = 2rj} and integers ki, ii, with ri — m — n < 
ki < £i < Ti + m, we have 

( q 

P j O {the line {z = 2r{\ has an x-gap D 



= det (l - [x[ fciA ](x i )K niP (2r i ,x i ;2r i ,x i )x[fc i ^](^)] 1 < iJ < g 
with a kernel K njP . 

The theorem was proved in [2] where the K-particles were called outlier 
particles. 

The dot particles of the L-process satisfy the following interlacing pattern. 

Proposition 1.3 For the "L-process, the lines £ = 2s contain blue dots and 
red dots, according to the following interlacing patterns, with varying num- 
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Figure 6: The red and blue L particles. The blue L particles correspond to 
Aztec diamond A while the red particles correspond to Aztec diamond B see 
Figure [T} The numbers represent the total number of particles on each line 
i = 2s 



hers: 



lines C, = 2s 


G Diamond 


# of blue and red dots 


< s < n — p 


e A 


n — s blue dots 


s = n — p 


e A 


p blue dots 


' n — s blue dots 


n — p < s < n 


e An B 


p dots with < 


to the right of 
s — n + p red dots 
k for each s 


s = n 


e An B 


p red dots 




n < s <2n — p 


e B 


p + s — n red dots 



with interlacing of the blue dots and interlacing of the red dots, with regard 
to the 7]- coordinate; in the overlap of the two diamonds, the dots interlace as 
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Figure 7: The interlacing system of blue and red dots. The p + 1 lines, 
{£ = 2s} <E AnB and {£ = 2(n — p)} 6 A, contain p dots. All the other lines 
contain more dots. See Proposition |1.3| for more details on the interlacing. 



well, with the right most dot on the line £ = 2s being to the right of the right 
most dot on the line £ = 2s + 2; also the left most dot on the line £ = 2s + 2 
is to the left of the left most dot on the line £ = 2s. Notice that the overlap 
contains p lines (through black squares) £ = 2s with n — p < s < n. 

Figure [7] shows the above proposition schematically for the example tiling 
in Figure [3j The proposition will be proved in section |4j 

It is interesting to notice that, passing from the L njP -process to the K np - 
process, the dot is maintained in the horizontal dominos, whereas a dot in an 
East domino gets replaced by a dot in a West domino; compare the pictures 
given in Figure [5] 

Recall the K niP -point process is a process of dots along the lines rj = 2k—l 
or what is the same z = 2k. For the sake of the main theorem below, the 
point of view will be switched around: namely, the K n p -process induces a 
determinantal process of dots along the lines £ = up to £ = 2(2n — p), 
inherited from the dots on the lines z = 2k. As mentioned, this process 
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can be obtained from the L n p -process by keeping the dots belonging to the 
horizontal domino's and moving the dots from the vertical domino's with 
a black square below to the vertical ones with a black square above; see 
Figure |5j 



1.3 The Tacnode GUE-minor kernel and the main 
Theorem 

We now define a new kernel, the coupled GUE-minor kernel, depending on 
two parameters /3, p: 

=K^{ Ul ,P- yi ; u 2 ,0-y 2 ) (2) 

+ 2/(1 - X?(\, «));i p A^(k), Bjf-^A)) 
x ' >-p 

where K minor (n, x\ n', x') is the GUE-minor kernel, defined for n, n' e Z, rather 

than N. This kernel will appear below as the appropriate scaling limit of the 

L-particle kernel. Here and below we shall define functions, which involve 

integration over small circles r and an imaginary line L := + + iM f; the 

line L needs to be always to the right of the contour IV Define 

K minor (n, x; ri, x') := - l n>n ,2 n - n 'U n - n ' \x - x') 

2 f r dw e - z2+2zx w n ' 
+ ~, — ^77 / dz 



(2vri) 2 7 ro J L w- z e - w2+2wx> z n ' 
with M. m (z) defined for m > 1 as 

„m— 1 

M m (z) := — 1, >0 . 

[m — ly. 

Kernel Q contains the functions: 

/ d( f du e-^+^K 
X (A, k) := - _ - g _ 2w2+4/3w ^ A+1 

* 1 J ' / ro (2vri) 2 7 L C-w e~^ 2 +^ J L 2vri C +K+1 



BJ-»(A) := 



r 



(2m) 2 J L (-co e-" 2 - 2 ^ o;- u J ro 2vri e u 2 -My+m ' 



3 The subscript > — p refers to the space £ 2 (—p, . . . , 00). 
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To be precise in the limit theorems we should replace l z > by 1 2>0 + §lz=o 



in the case of the L-process, (1.4) below, and by l z>0 in the case of the 
K-process, (1.5) below. Since these changes do not affect the limiting point 
process we will ignore this fine point. Some properties of the kernel are given 
in section [3] Notice that the scaling in the Theorems below could have been 
derived from the scaling used in the limit of the IK-kernel to the tacnode 
process, combined with the way the weight a — > 1. 
The main statement of the paper reads as follows. 

Theorem 1.4 Let the size n = 2t + e, e G {0,1}, of the diamonds go to 
infinity, while keeping the overlap p = n — 2m finite and, and letting the 
weight of the vertical domino 's a — > 1 as 

with G R fixed. 



The coordinates (£, r/) are scaled as follows, 

& = At + 2e - 2u h r]i = 2t + 2[yiVi] - 1, with Ui G Z , y { G R. 
With this scaling, the following limit holds: 

hm (- a )^-^)/ 2 (-Vt)« 1 -^)/ 2 L ni p(ei^i;e2,^)V^ 

n— >co 

= K^,(ui,yi;u2,Jte)- 

We also have that the rescaled h-particle process converges weakly to the 
determinantal point process given by the tacnode GUE-minor kernel. 



Remember the remark at the end of subsection L2) The K-process in 



duces a process of particles along the consecutive lines {£ = 2s}. This is to 
say the roles of 2r and x in the kernel K n p get reversed from the point of 
view of scaling: the variable 2r turns into the discrete variable u and the 
variables x into the continuous variable y. The simulations of Figure [9] show 
that the lines {£ = 2s} belonging to the overlap A D B contain long dense 
stretches of K-particles. But performing a random thinning, one nevertheless 
is led in the limit to a point process kernel, which turns out to be the same 
tacnode GUE-minor kernel, except for some shift. This is the content of the 
next Theorem. 
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Theorem 1.5 Let a, n and p be as in the previous theorem and consider the 
same scaling as above, but expressed in the (x, z)- coordinates, using the map 

2 Xi = -p + 2 Ui -e + 2[ yi Vi}, 2r t = 2t + 2[ yi Vt}. 

Let the K-particles be thinned out at the rate p n — 1 — 2/\/t. We then have 
the following limit: 

lim (1 -p n )a r2 " ri (v / t) a;i - :C2 - ri+r2 (-l) :El - X2 K niP (2r 1 ,x 1 ;2r 2 ,x 2 )v / t 

n— ¥oo 

= K% p (u 2 + l,y 2 -u 1 ,y 1 ). 

Interpreted as a weak limit of a point process this means that if we thin the 
K-process by removing eachWL-particle independently with probability p n , then 
the resulting point process converges weakly to a determinantal point process 



given by the correlation kernel on the right hand side of (1.5). 



Notice the kernel K^ ac p is the same as the one in Theorem 1.4, except for 
the shift in u 2 and the flip u x -h- u 2 and y\ <->■ y 2 . 

The geometry of the level curves for the double Aztec diamond, when 
n — > oo, looks as in Figure [8] below. The K and L particle processes have also 
been plotted for this simulation and is found in Figure [9] below. As n — > oo, 
the particles take continuous values on each line; they are constrained by the 



same interlacing as in Proposition 1.3 



2 The kernel for the L-process, via Kasteleyn 
2.1 The L-particle process 



The kernel L np for the L-process is given by formula (1.1), i.e. 



L» >P (£i, w & v*) =(! + « 2 )4 0) (6, m\ 6, m) 

- (1 + a 2 )(((J - K n );l p+ A Cum )(k), B^ m (k))^. P+1 , 
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Figure 8: A simulation of a double Aztec diamond with n = 100 and p = 4 
with the weight a of vertical and horizontal tiles equal to 1. Both figures are 
rotated by 7r/4 counter-clockwise. For this simulation, the top figure shows 
the underlying domino tiling while the bottom figure shows the level lines. 
The simulation was made using the generalized domino shuffle |23j. 
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Figure 9: The top picture shows the blue and red L-particles for the simu- 
lation of Figure [8j a particle is created each time a line {£ = 2s} traverses a 
blue or green domino, as in Figure [8] (or, alternatively, crosses a level line); 
then blue particles belong to an A-level line and a red particles to a B-level 
line. The bottom picture shows the K-particles for the same simulation in 
Figure [8j a particle is created each time a line {r/ = 2r + 1} traverses a red 
or green domino, as in Figure M (or, alternatively, crosses a level line). 
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with 



+ 



f dw (l + aw)^ 1 -^^-! tl _ (a 
" (6<£a) 7r a 27i ( w -a)to~^)/2+i w 

/" rfw (l + az)(" 1 - 1 )/ 2 (z-a) n -( r ' 1+1 )/ 2 w n -« 2 / 2 



r 0l , 



ia (2?ri) 2 y roa2 w-z(l + a^)te +1 )/ 2 (w - a)"-(^-i)/2 z n-6/2 
x ) i+fc I I dz z n - j (l + aw) n (w - a) n+1 



h ,_{-V J+ f J" dz z n -\l + aw) n {w - a) n+L 
Ah ] " ~W¥ Jv , a ^ w^il + azHz-a)^ 



(27ri) 2 J ro a 7r ,a,* ^-z (1 + aw) n (w - a) n+1 z n -&/ 2 

_\\k j ^ei/2-fc 



2vri J roo (1 + - a)^+ 3 )/ 2 ^ 

(-l) fc /" /" cfe (1 + ato) n (w - a) n+1 z n ~« 2 / 2 

«) - ^ rfw — (l - az y m+1)/2 ( z _ a ) n - te _i)/2 w „+i-fe 

+ hil! / (1 + az)™-^ +1 )/ 2 (^ - a )te+i)/2^-i-6/2 ^ 
2?ri /r ,„ 

(4) 

As was shown in [2] the IK-particles form a determinantal point process 
and the kernel K n p for the IK-process is given in (z, x) coordinates by 

(-lf-^K^r, x- 2s, y) = K° (2r, x; 2s, y) 
where K n (j, k) is defined as before in Q and where 

K° (2r, x; 2s, y) = K°;f zt (2(n - r + 1), m - x + I; 2{n - s + 1), m - y + l) 
= -ls<r.(-l) x "^2r, 2s (a;, y) + S(2r, x; 2s, y) 

f_yk-x r r dv v x+m ( 1 + ay y( 1 _ ^n-s+l-e! 



>n-p+l 



ax,28+ ci (k) ■= /0 .x 9 f du 



(2vri) 2 J Toa J Toau u-v (i + «)»(l-a)«+i 

-l) fc - x / / cfu ^ +m (l + at») s (l- «)»»— 

aw A 



(27ri) 2 J ro a 7 r0iBj . « - « (l + aii )»(l_ft)n+l 

dv- 



2m J v (l + tM ;)«-(l -*)•+« 
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y ' 2r+€A) • (2vri)2 f Voa dU f Voau v - u uv^ (1 + any {I - ft)»-H-i-* 
(-l) 1 -" jf jf v 1 (1 + au) n (l - ^ n+1 



i 2 ^) 2 Jr , a Jr , a , v v-u v,y+™+ 1 (1 + auY {I 
i) l -y f , (1 + av) n - r (l - 2 ) r+£2 



\n— r+l— e2 



27Ti /p, ^jz+m-I+l 

1 0,o 



f y f f dv 

S(2r + ei, x; 2s + e 2 , y) := v , n & du A 



(27ri) 2 7p n J T v-u 

\ / ^ 1 O.a ^ 1 O.a.u 



u y-m n + av )r(l _ a)n-r+l-ei 



r0i- 2™ (1 - 

(5) 

A single Aztec diamond of size n leads to a determinantal process as well 
(see [15]), for which the kernel is given by the following expression: 

K° aeAzt {2r,x;2s,y) 

(-iy-y r r dv v x (l + au) n ~ s (i - f y / 

(2tti) 2 J 7r3 J 7r2 d-h 1 f (1 + av) n r {l - -) r 

It is not immediate to go from knowing the kernel for the K-particle 
process to the kernel for the L-particles process. To do so we will use the 
fact that we can get the inverse Kasteleyn matrix for the dimer version of 
Double Aztec diamond. The inverse Kasteleyn matrix will be explained in 
terms of the kernel ]K n p . Using the inverse Kasteleyn matrix it is possible to 
show that the L-particles form a determinantal point process and compute 
the kernel. 



2.2 The Kasteleyn Matrix 

Suppose that G = (V, E) is a bipartite graph. A dimer is an edge and a 
dimer covering is a subset of edges such that each vertex is incident to only 
one edge. The dual of the double Aztec diamond is a subset of the square 
grid graph with a certain boundary condition while a domino tiling of the 
double Aztec diamond is a dimer covering of its dual graph. Kasteleyn, in 
|19j . introduced a matrix, later named the Kasteleyn matrix which one can 
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use to compute the number of domino tilings of the graph. Since the graph 
in this paper is bipartite, the Kasteleyn matrix is a type of signed weighted 
(possibly complex entries) adjacency matrix with rows indexed by the black 
vertices and columns indexed by the white vertices. The sign of the entries 
is chosen so that the product of the entries of the Kasteleyn matrix for edges 
surrounding each face is negative. This is called the Kasteleyn orientation. 
We will describe the Kasteleyn matrix for the double Aztec diamond below 
but first we state Kasteleyn's theorem for bipartite graphs and Kenyon's 
formula [20j. 

Suppose that K denotes the Kasteleyn matrix for a finite bipartite graph 

G. 

Theorem 2.1 (||19|) The number of weighted dimer coverings of G is equal 
to |detiT| 

Suppose that E = {ej}" =1 are a collection of distinct edges with e« = (bi,Wi), 
where 6j and tOj denote black and white vertices. 

Theorem 2.2 ( [20J ) The dimer s form a determinantal point process on the 
edges of G with correlation kernel given by 

L(ei,ej) = K(b h w i )K' 1 (w h b j ) 

where K(b,w) = K bw and K~" 1 (w,b) = (K _1 ) w6 

The above theorem means that by knowing the inverse of the Kasteleyn 
matrix, which we call the inverse Kasteleyn matrix, we can derive the corre- 
lation kernel of the dominos. We can now introduce the Kasteleyn matrix of 
the double Aztec diamond. 

Let 

W = {(x 1 ,x 2 ) : x x e 2Z + 1, x 2 £ 2Z} 
denote the set of white vertices and let 

B = {( Xl , x 2 ) : Xx e 2Z, x 2 G 2Z + 1} 

denote the set of black vertices. The dual graph of the double Aztec diamond 
written in (£, if) co-ordinates has white vertices given by 

// l<x 1 <2(2m + 7i) + l,0<z 2 <2(7i-l) 1 

Wad = < {xx,x 2 ) e W : n ' ' > 

[ or 1 < xx < 2n — 1, x 2 = 2n j 
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(0,3) 1 (1,4) 
-O 



2( n +2m) + 1 




(1,0) -1 (2,1) 



Figure 10: The left hand figure shows the dual graph of the double Aztec 
diamond with n = 8 and m = 4 in the (£, rj) co-ordinates. The right hand 
figure shows the weights, Kasteleyn orientation, black and white vertices for 
the two most left squares, with (£, ^-coordinates. 



and has black vertices given by 

< xi < 2(2m + n), 1 < x 2 < In - 1 



ad ^{x 1 ,x 2 )eB: or 2 (2m + 1) < xi < 2(2m + n),x 2 = -1 
We denote by A n ^ m to be vertex set of the dual graph of the double Aztec 



diamond with (£, rj) co-ordinates. Figure 10 shows the dual graph of the 
double Aztec diamond with n = 8 and m — 4. 

Let K a denote the Kasteleyn matrix for the double Aztec diamond with 
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( (_l)-(fri+6a+l)/2 a ( r ) if W = b + e r 

AT a (b,w)=^ (_i)-(&i+fe-i)/2 a ( r ) if w = b-e r (6) 
[ otherwise 

with d = (1,1), e 2 = (-1,1), b = (bx,b 2 ) G B, r = 1,2, a(l) = 1 and 
a (2) = — ai. The choice of sign for the entries of the matrix is the same 
as [3]. These are chosen so that entries of the inverse Kasteleyn matrix are 
discrete analytic functions when a = 1. 

Theorem 2.3 The entries of inverse Kasteleyn matrix for the double Aztec 
diamond, K a , defined by ^ are given by 



K 1 

= -(-!)(— ^- b ^K n , p fb 2 + 1, b2 ~ bl+ 2 2m + 1 ;w 2 + 1, 



b 2 — b\ + 2m + 1 1 w 2 — Wi + 2m + 1 

2 

(7) 



where as before n — p = 2m. 

In other words, using ([TJ, (x,z) -H- (^,77), 

K\v,b) = -(-1)(— 2+6l - fe+2 )/ 4 K n , p (^(b),x(b)^(w),x(b)) 

so essentially the two kernels are the transpose of each other modulo the 
co-ordinate transformation ([!]). 

The proof of this theorem involves the characterization of K~ x : Ka.K^ 1 = 
I and is given in Section [6] In order to find such a formula for K~ x we used 
a guess following the approach in j3]. More explicitly, using the IK particle 
correlation kernel one can compute the joint probabilities of K-particles. As 
these particles correspond to east and south dominos, this joint probability 
should be equal to a corresponding formula written in terms of the inverse 
Kasteleyn matrix by using Theorem |2.2| These two sides can be compared 
which gives a guess for the inverse Kasteleyn matrix in terms of the IK particle 
correlation kernel and leads to the formula in the theorem. 



Since we now have the inverse Kasteleyn matrix we can use Theorem 2.2 



to prove theorem 1.1 The basic observation is that we have an L-particle 



at a black vertex b if and only if a dimer covers the edge (6, b + e±) or the 



edge (b,b — e 2 ). By using theorem 2.3 we can compute the probability of 
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seeing L-particles at given black vertices b\,...,b£ by summing over all the 
possibilities for the dimers and using theorem |2.2| and thus deduce the L 
kernel, see section [7] for the details. 

For general weights and boundary conditions of the square grid, if we 
define a point process on the black vertices such that a particle is present 
at a black vertex iff a dimer is incident to that black vertex, from Theorem 



2.2 , we can recover the particle correlation kernel provided we know the 
inverse Kasteleyn matrix of the model. In general, the reverse, i.e. to go 
from the particle correlation kernel to the inverse Kasteleyn matrix, is quite 
complicated. However, for the double Aztec diamond, we were able to express 
the inverse Kasteleyn matrix in terms of the K-kernel. There should also be 
an analogous formula for the inverse Kasteleyn matrix in terms of the L- 
kernel. This formula could then be used to give the K-kernel in terms of the 
L-kernel. 



3 The tacnode GUE-Minor kernel and its 
symmetry 

Recall from ^ the tacnode GUE-Minor kernel, about which we show the 
following. 

Proposition 3.1 The kernel K.p c (ui, y%\ u 2 , y 2 ) is invariant under the invo- 
lution 

ui -H- p-u 2 , yx -H- -y 2 , 
and is a finite rank perturbation of the GUE-minor kernel, as follows: 

(J K mmor ( Ml ,/3 - yi - u 2 ,(3-y 2 ) 

max(p— l,p— 1— U2) 

+ 2 ((l-^(A-p,K-p))- 1 <^)(A-p)^f^(A-p) 

A=0 

= ) K mmor (p -u 2 ,P + y 2] p- Ul ,P + yi) 

max(p— — 1) 

+ 2 J2 ((! - - Pi K ~ PT'K-u^) (A - p)S^r'(A - P) 

A=0 
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Remark: This symmetry (3.1) is not surprising, since it corresponds to the 



symmetry of the geometry of the double Aztec diamond. 



Proof: In order to prove this statement we need the following functions, 
G(X) = [ —e 2 " 2 - 4 ^u- x - 2 , g y (k) - f ^""-^-vK 



, e ^ tu 

l 2vr« J L 2m 

1 ; J To 2vr<- vK } J Fo 2ttzC- a 

and the corresponding operators 

9(«, a)f(a) := ^ G ( K ~ P + a )f( a ) 
J£(A, a)f(a) := F(A - p + a )f( a ) 

a>0 



The kernel % /B and its resolvent, and the functions A^ y (n) and !Bf ,y (A) as in 
^ can then be expressed as follows, 

X p (X, k) = G ( X + a ) H ( a + K )> X p( A - K ) %0 ( X ~ P, K ~ P), 

a>0 

x p = sh , % T P = rcs, with g T = g , w T = jc, 

1 + X(X, «) := (1 - X p (A, /t))" 1 = X a p 

a>0 

- p) = ff-itt (« - P + mi) - ^ G(k - p + a)/^ (a - u x ) 

= -P + u l) - 9(K, - «l) 

^f^A - P) = /i-, 2 (A - p + « 2 ) - ]T H(X - p + a)^ 2 (a - u 2 ) 
= h- y2 (X -p + u 2 )- IK(A, •)&&(■ - "2) 
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>0 



>0) 



>0 



>0 



Using the definition ([2]) of the kernel and the expressions (J3]), we have the 
following identities: 

= -1^2^-^-^ (y 2 - y x ) + Y,9yA<* - u 2 )h yi {a - «i) 

a>0 

+ ^(1 + K))g- Vl (K - p + ui) , h- y2 (X ~ P + U v) 
+ a)g y2 (a - u 2 ) , (1 + ^(A, k))S(k, a)h yi (a - «i) 

- ^(1 + K(A, K))g_ yi (n- p + ux) , :K(A,a)^ 2 (a -u 2 )\ 

- ^(1 + CR(A, «))S(«, - u x ) , h_ y2 (\ - p + 1*2)^ 

= -i Ml>U2 2" i -" 2 - 1 e« i -" 2 (y 2 -yi) 
+ (g- y A K - p + u i) , (1 + ft T (A, K))h_ y2 (K - p + w 2 )y 

+ (^ 2 (a- M2 ) , (i + ^ T (i + ^)g)/ ij/1 (a- Ml )y 

- + tyg- yi (K - p + ui) , ^ 2 (a - M 2)) >o 

- ((1 + %)$h yi {{a - Mi) , h^ y2 (\ - p + W2)) >Q . 

Given the involution (3.1), all terms in the last expression (|3]) are self-dual, 
except that the second and third terms interchange, because of the operator 
identity (see (|3|) 

1 + JC T (1 + R)S = 1 + 5C(1 + 31) S = 1 + £ T 

and the self-adjointness of CK(1 + "Jl) and (1 + 3?) 9- 

To prove the second statement ( 3.1[ ) on finite perturbation, one notices 
from ([3]) that the double integral in CBf ,y (A) equals for A > 0, since the 
integrand as a function of ( has no pole at and, similarly the single integral 
equals for u + A > 0. Thus S^ 2_/3 (A — p) = for A > p and for A > p — u 2 . 

This proves from pj), t he first equality =. The second equality = is obtained 
by the involution (3.1). ■ 



>o 
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4 Interlacing pattern of the L-process 

In this section we prove the interlacing properties as explained in Proposition 



1.3 We first need the following Lemma (remember p 



n 



2m). 



Lemma 4.1 The total number of dots along the line £ = 2i equals the dif- 
ference of height Ah between the extreme points of that line, as is given by 
the boundary values of the height function: [j 



lines £ = 2i 


hleft 


h 

' v r~ight 


Ah 


# of blue and red dots 


< i < n — p 


n 


i 


n — i 


n — i blue dots 












p + i — n red dots 


n — p < i < n 


p + i 


i 


P 


1 


to the left of 












n — i blue dots 


n < i <2n — p 


p + i 


n 


p+i — n 


p + i — n red dots 



Proof: The statement on the first row in the table above follows from the 
fact that the height h along the lines £ = 2i for < i < n — p decreases 
from h lctt = n to h right = i for < i < n — p (going from left to right) and 
from the fact that each decrease of height by 1 produces a dot. The same 
statement holds for the range on the third line of the table by the obvious 
symmetry consisting of flipping the figure about the middle of the £-axis and 
the middle of the //-axis. Also note that the heights of the S-level curves 
range over the half-integers from n + 1/2 to 2n — 1/2. Therefore the lines 
£ = 2i for < i < n — p, which have height at most n, will never intersect 
those lower- level lines and vice- versa, showing that on the first line (resp. 
last line) of the table above only blue (red resp.) dots appear. 

In the overlap region of the two diamonds, the boundary values of the 
height function show that h left , h Iight and Ah = h left — h Iight is as indicated 
in the table. Moreover, since the height of the 5-level curves is > n + 1/2 
and the height of the A-level curves is < n — 1/2, the red dots are all to the 
left of the blue dots along the lines £ = 2{n — p) up to 2n, with numbers as 
indicated in the table. 

Lemma 4. 2 Let the lines £ = 2k and £ = 2k + 2 for < k < n — 1 have £—1 
dots starting from the right boundary. Then the £ th dot on the line £ = 2k + 2 
must be to the left of or coincide with the i th dot on the line £ = 2k. 

a dot-particle x is to the right of a dot-particle y means rj(x) > ij(y). 
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Proof: Note that the right most point of the double Aztec diamond on the 
line £ = 2i has height i, provided < i < n. Therefore, if there are £ dots 
on the line £ = 2k + 2, counting from the right hand boundary, then the left- 
lower vertex of the square A', containing the £th dot, has height k + £ + 1; 



see Figure [TTJ Consider two cases: 

(i) assuming no dot in the corresponding square A on the line £ = 2k, 
then the only way to cover the squares A and A' with domino's such that 
A' carries a dot and not A, is given by the four upper configurations of 



Figure 12; putting in the heights forces the height of the left-lower and right- 



upper vertices of the square A to be k + £ as indicated in Figure 11 This 
shows there must be £ dots on the line £ = 2k strictly to the right of the 
square A. So, the £ th dot on the line £ = 2k + 2 must strictly be to the left 
of the £ th dot on the line £ = 2k, at least if A contains no dot. 

(ii) Assume a dot in A on the line £ = 2k; the only way for this to occur is 
given by the four lower configurations of Figure [12} From them one deduces 
that, if the height of the lower- left corner of A' is k + £ + 1, then the height 
of the lower- left corner of A must be k + £ or k + £ + 1. In the former case 
(i.e., k + £), the dots in A and A' are the £th ones from the right, proving 
the claim; in the latter case (i.e., k + £ + 1), the dot in A is the £ + 1st one 
and the dot in A' the £th one. So, the £th dot on the line £ = 2k is to the 
right of the £th one on the line £ = 2k + 2. ■ 



Proof of Proposition L3 Consider two consecutive lines £ = 2a and £ 



2a + 2 through blue dots, with the squares A and B, containing each a dot 



and no dot in between A and B; see Figure [13} This is to say, the level of the 
line £ = 2a goes down from k to k — 1 within square A, stays flat in between 
A and B and then goes down from k — 1 to k — 2 within square B. We now 
consider the line £ = 2a + 2 between the two corresponding squares A' and 

with same rj coordinates as A and B respectively. 
We show there must be at least one dot in between the squares A' and B' , 
possibly including A 1 or B'. One checks there are exactly six configurations 



with a dot in the upper-left square; see Figure [14} Superimposing any of 
the four upper configurations on (A, A') or (B,B') will give a dot in A' or 
B' . Assuming no dot, neither at A 1 , nor at B' , the configuration (A, A') or 



(B, B') at Figure 13 can be covered by any combination of configurations (I) 
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k 



k. + l 




Figure 11: Assume the £th dot on the line £ = 2k + 2 (counted from the 
right) appears in A', and assume no dot in the square A, then the height of 
A must be as indicated. 



i + 1 



i + 1 



i + 1 



i - 1 



/' + 1 











■ 




5 


i - 


- i 



i + 1 



Figure 12: The four upper configurations are the only coverings of A and A' 
of Figure 11, with A' carrying a dot and not A. The four lower configurations 
are the only coverings of A and A', with both A and A' carrying a dot. 
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fc-2 



D 



• 




- 1 





B' 



A 






k-i 







U' 



Figure 13: Between the two gray squares labeled A and B the height function 
stays constant; therefore the line between A and B contains no dots. 







i i i 



i i +i 



/' + 1 i + 1 /' + 1 



' - ! - I ! \ ' M J + 1 II 

7+1 7+1 7+1 



Figure 14: Let the squares A and A' (as in Figure 13 ) each contain a dot, then 



the four upper figures are the only possible covers of (A, A'). If A contains a 
dot and A' does not, then the two lower figures are the only possible coverings. 
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and (II) in Figure 14 Indeed 



{A, A') 
I(i=k-1) 
I(i=k-1) 
II(j=k-l) 
II(j=k-l) 



(B,B>) 
I(i=k-2) 
II(j=k-2) 
I(i=k-2) 
II(j=k-2) 



In all four cases, the difference in height between the lower- left vertex of A', 
having height k, and the upper-right vertex of B' , having height k — 1, will 
always = 1, thus creating a jump in between and thus one dot. So in all 
cases, there will be at least one dot in one of the squares on the segment 
(A',B'), including possibly on the extremities. 



Finally, this fact together with Lemma 4.1 on the number of blue and red 



dots and Lemma 4.2 imply the interlacing, with regard to the ^-coordinate. 



5 Scaling limit of the L and K -processes 



In this section we will prove theorem 1.4 and theorem 1.5 Let 5 G {0,1}, 
where 5 = will correspond to the L-kernel and 5 = 1 to the K-kernel. We 
will ignore the integer parts in the scaling (1.4) and (1.5). This makes no 
essential difference but simplifies the notation. Set 



Then the prefactor in (1.4) for the L-kernel can be written 



and the prefactor for the K-kernel in (1.5) is 



,r2— ri 



xi—X2+r2-r± i 



\X1-X2 



(Vt) 

2-J/l)v / t^_ a ^(2/l-«/2)v / t^_ v /^^«l-(«2 + l)^_ v /^^ 



,2(2/: 



,2{y 2 -yi)Vt 



ft,i{u 2 + l,y 2 ; ui,yi). 
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Define 

C 2t+e, P ,5{ u ^y^ u ^V^) = -ft,s( u i>yh u 2l y 2 )(l + a 2 ) ((1 - 5)l U2<Ul +5t yi<y2 ) 



x 



2vri 



To 



(1 + aw) {yi - y2)Vi - 1+s 

a (w - a)^-y^ t+1 - & ' 



U2—U1 



dw, 



^2t+e,p,6 



X 



(u 1 ,y 1 ;u 2 ,y 2 ) = /t,5(wi,yi;w 2 ,y 2 ) 
( 1 + fl2 ) / f dw w U2 (l + az) t+yiVi - 1+s (z - ay-y iV ~ t+e+s 



(27Ti) 



ro,a 



r ,a, z w - 2 z"i(l + awY+y^iw - a)*-^ 2v/ * +1+e 



and 



,yi;u 2 ,y 2 ) 

= — /*,* (til, 2/15 1*2, 2/2) (1 + ^ 



X M(J -^2t+e)>2t+ e -p+l y 44t+2e-2«i,2t+2yiv / t-l,<5)0') 



5 



4t+2e-2M 2 ,2t+2y 2 v / i- 



>2f+e-p+l 



where 

^4t+2e-2ni,2t+2j/iv / t-l,<5(^) 
(-1)' 



+ 



2vri 



2t+e— Ui — k 



dw 



(27ri) 



dz 



T 0az w - Z z u ^(l + aw) 2t+t (w - a) 2t+e+l 

which is a slight modification of A^ ltVl (k) in Q, and where B^ 2>m {k) and 
K 2t + t are as given in Q. With these definitions it follows from (2.1) that 

(-a)to-^)/ 2 (-^)^- 6)/2 L 2t+e> p(ei^i;e 2 ^2)v^ 



(1) 



2/ + . /mi T~ C2t+e,p,0 

,yi,u 2 ,y 2 ), 



»(3) 



if we have the scaling (1.4). Similarly, it follows from (2.1) that 



\x\-X2+ri-r\ , 



1 + a 



-l) x ^K 2t+e>p (2r u x 1 ;2r 2 ,x 2 ) 
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Figure 15: The contour paths 



We will now use ([5]) to prove (1.4). The proof of (1.5) from QSJ) is completely 
analogous since the change from 5 = to 5 = 1 has no effect in the limit. 
Note that a 2{yi ~ y2)Vi e my2 ~ yi) as t -»■ oo since a = 1 - (3/y/t. We see that 
(1.4) follows from 

3 

SE C 2?+v,o( M i'2/ii "2,2/2) = e 2 ^ yi ~ y2) K 0tP (ui,yi] u 2 ,y 2 ). 



i=i 



Let Ci be the positively oriented unit circle and let C 2 = C' 2 + ^2> where 
C 2 consists of two infinite line segments t — > j3 + it, t e (— oo, —2] U [2, oo), 
and C 2 is a smooth curve that goes from (3 — 2i to (3 + 2i to the right of Ci, 



see Figure 15 



Let C 2 be C 2 reflected through the origin. 
Set 



and 



*Im(C) 
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so that F X)t (C) = F 0tt (C)G Xtt (C). Also, we write 

9x A0 = e^-<\ and //3 (C) = e 2 ^ 2 . 
The next lemma contains the estimates we need. 

Lemma 5.1 Fix A > 0, (3 e K and > 1. There is a to cmd a constant C 
such that for all t > to, x e [—A, A], s G K and C G Ci U 63 we have the 
following estimates 

1 1 

< 



\F x , t (P + is)\ ~ l + s™/2 k kV 



and 



G x ,t(0 1 

9xA0 

FoAO 1 



MO 



< 



< 



c_ 
c 



Proof: We have, since a — 1 — (3/ that 

t+xVi 

\F x , t (P + \s)\ 



a- 1 - 



/3 + is 



Vi 



Now, 



and 



13 + is 



Vi 



1 + 



/3 + is 



Vt 



t—x\fi 



.1 + 



Vi 



( ^7r-^V + T> 1 + 



when t is large enough. Thus, 

/ 2\ §(*+Vt) f J\ |('-^\/*) / „2\* 

|J^(/* + i*)|>(l + T J (l+ T j =(l + 

By the binomial theorem, 



t 



fc=i 



> 1 + 



(*/2)« 



,2fc 



1 + 



^2fc 



2 fc A;! 
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if k < t/2. This proves (5.1). Note that Ci U C 2 ' is a fixed compact set. The 



estimates (5.1) and (5.1) follow from the inequalities 



Mlog(a + C/Vi) - x(( - 0) 



< 



C 

7r 



X 



s/ilogia' 1 - CM) - x{fi - C) 



< 



C 



and 



tlog^a- 1 - CM) (a + CM) " (2C/5 - C 2 ) 



< 



for sufficiently large t, which in turn follow from Taylor's theorem. 



Consider first C 



(i) 

2t+e,pfl- 



The case y\ = y 2 is special. In this case we obtain 



c S+ eiP ,o( w i>?/i; w 2>ya) 



-(-V*) tla ~ ,tl V'*l 



U2<Ul 



-Vt 



U2<Ul ' 



1 + a 2 
2vri 



a + 1/a 
2vri 



it;' 



(l + aw)(w-a) 



_ x (w + l/a)(w - a) 



dw 



-(-Vt) u *- u 'Vil U2<Ul (-l/a) 



U2 — Ml 



In the second inequality we deformed the contour through infinity to a con- 
tour surrounding —1/a. If «2 = u\ — 1 this equals —a which goes to —1 as 
£ — 7- 00. If k, 2 < U\ — 1 the last expression goes to as t — > 00. 

If yi > y 2 then deforming the contour to r_iy shows that C^]^ „ = 0. 
Assume that y 1 < y 2 . Then, for large enough t, since w = a is not a pole, 



C'i+e.p.oK^i;^,^) 



-1 



U2<Ul ' 



L U2<Ul 



1 + a 2 
2m 

1 + a 2 
2m 



r 



a — w 



l-wy/i) 



U2—U\ 



(1 + aw)[a — u;) 



«2— Wi 



(1 -au}/y/i)(a + u/y/i) 



dio 
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by the change of variables w = —u/y/i. It now follows from lemma 
lim C£l (u 1 ,y l ;u 2 ,y 2 ) 



5.1 



that 



_H —I Jto-mKv-ft duJ 



U2<Ul 2niJ ei ~ 



U2 — 1 



I „2/3(y 1 -y2)nu 1 -u 2 (y^__y}y 

(Mi - U 2 - 1)! 



Thus, for all y±,y 2 , 



lim (u 1 ,y 1 ;u2,y 2 ) 



t— >oo 



\Ml-«2-l 



— U-u 1 =u 2 +l JL yi=y2 JL u 2 <u 1 ^y 1 <y 2 K * _ u _ 

= _ e 2p{y 1 -y 2 ) 2 u 1 -u 2H u 1 -u 2 ^ y2 _ y ^ 

where 

ytn— 1 i 

for m > 1. 

Consider now C^t+e _o- We can write, using the fact the z contour has 
no pole at z = a for t large, and by completing the G 2 /\/i contour with an 
infinite semi-circle on the right, 



c£ ) +e,p,o(u 1 ,y 1 ;u2,y 2 ) = jt^^ <£ dz(-Vi) 

dw (-wVi) U2 (a- 1 + z) t+yiy/i (a - z y- VlVi+e 



e2~M w ~ z (-zy/i) Ul (a -1 + w) t+ y^(a - w )*-2/2^+e (l + a z)(a - w) 
1 + a2 I d( I duj " U2F viA0 (a + C/Viy 



(27ri)2/ ei Ve 2 ^-CC Ml ^ 2 ,tM(l-aCM)(a + a;M) 1+s ' 
It now follows from lemma IB~T1 that 

lim C£l (u!,yi;u 2 ,y 2 ) 



t—>-oo 



du u^f p {0g Vl A0 



(2vri) 2 / Cl Je 2 uj-( C^Z/jMswH 

2T r din w u 2 p -C 2 +2(/9-yi)C 



(2vri) 2 / ei 7e 2 ^ - C C Ul e~^+ 2 ^-y^ ' 
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We choose k in (5.1) so that 2k > u 2 , which gives a uniform t-independent 
upper bound on C 2 . Combining ^ and ^ we see that 

lim (Co + ^^oK"!. Vi; «a, 2/2) = e 2 ^ 1 -« a )K minOT ( Ul , /3-y i; « 2) P~y 2 ). 

Next, consider 



00 

2j S 4t+2e-2ii 2 ,2t+2 2 , 2 v / I-l(^ + 2t + e - p + 1) 



00 

X 

K,\=0 

x (7 - K 2t+ e)>o(A + 2t + e- p + l,re + 2* + e- p + l) 

X y 44t+2e-2ui,2t+2^i V / t-l,o( K + 2t + e — p + 1), 



where A^ uVl>0 is given by (j5]), 

^4t+2e~2u 2 ,2t+2y 2 VI- 1 U ) 

i-iy 



2ixi 



(1 + a^)'-^ +e (^ - a ) t+y2Vi z U2 - 1 - 2t ' e+j dz 



It. a 



f dw f dz z U2 (l + aw) 2t+e (w - a) 2t+e 



(27ri) 2 J Foa J Toaw w-z «,«+H+i(l + az)*+w^(z - a)*-^+i+e 

and 

/ n (-iy +k f r f dz z 2t+ '-Hl + aw) 2t+ Hw-a) 2t ^ +l 
K 2t+t (j,k) = / dw 



Set 



(2yri) 2 J Toa J Foaw z-w w v+c~k+i(i + a^) 2 *+ e (z - a) 2 *+ 



X ^4t+2e-2M 1 ,2t+2j/ 1 ^-l,o( K + 2t + C — p + 1) , 



e+1 ' 



X B 4t+2e-2u 2 ,2t+2y 2 Vi-l( X + 2t + C - p + 1) 

and 

iW(A, «) = (Vt) K ~ x K 2t+ e(X + 2t + e-p+l,K + 2t + e-p + l). 
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Note that the matrix with elements 

(Vt) K -\l - K 2t+e )>l(\ + 2t + e-p+l,K + 2t + e-p+l) 
is the inverse of the matrix with elements 5 Kt \ — K 2t+e (X, K ). Thus, 

oo 

c £+ e ,pfi(ui,yi; u 2 ,y 2 ) = -(1 + a 2 ) B„(\)(I - K 2t+e )>lA Uuyu0 (K) 

re,A=0 

and we want to take the limit of this sum. (Note that the sum is finite even 
in the limit.) 

Rewriting in the same way as above we see from ^ that 



-re— ui+p— 1 



dcj 



27ri/ e2 F_ w , t (w) (1 - au/Viy+^a + u/Vi) 



d( 



du C- Ul F yi , t (0 



c/Viy 



(2vri) 2 / Cl Ve 2 oj - ( c^+^Fo^) 2 (1 - <M)(a + u/Vt) 1+e (l - aju/Jiy 



Using lemma 5.1 we can take the limit t — > oo to get 



1 



lim Aii,itt,o(«) = ,> ■ r t , ^ 

2tti / ea f p {u)g_ yu p{u) 



— re— wi+p— 1 



do; 



+ 



(27Ti) 



e 2 



Similarly we get 

B U 2,V2 (^) 



+ 



1 



2?ri 



da; 



— y2,* 



dC a; A -"F , t (a;) 2 (a + a;/^) 1+e (l - aa;/ v / t) e 



(2vri) 2 JT ei / e2 a; - C C^^CO (a + CM) 1+t 

and again by lemma 5A_ we find 

1 

^?jo.?;o( A) — ; 

t— >oo 



lim B U2 ,y 2 {\) . 

Z7T1 



U2+X—p 



e-i 



+ 



dC ^"fp^f 



(27ri) 2 / ei J C2 a;-CC- U2 MC)<w(C) 
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Finally, 

Sr ,* x 1 I j I d( u K ~ p Fo,t(*r (l-au/Viy{a + uj/Viy +1 
and we see from lemma 15.11 that 



(2vri)2/ ei / e2 C-^C A - p+1 ^(C) 2 
= X fi (X-p,K-p). 

It now follows from g, Q, Q and @ that 

Together with ([5]) this proves (J5|). 

It is not difficult to get t-independent bounds on the L-kernel using the 
same arguments as above and in this way we can show, in a standard man- 
ner, that the appropriate Fredholm determinant converges and obtain weak 
convergence of the L-particle point process. We will not enter into the details. 



6 Proof of the inverse Kasteleyn formula 

In this section we prove Theorem 2.3| We will use the fact that 



b 2 - b\ + 2m + 1 ^ w 2 - Wi + 2m + 1 
t &2 + 1, ; w 2 + 1, 

r^inUer fr , , b 2 - h + 2m + 1 W 2 - Wi + 2m + 1 
= ~KT ( 6 2 + 1, g ' W2 + ' 2 

where the kernel K 1 ^ 1 is the inlier kernel from [2], dual to lK njP . We will 
use the form of the inlier kernel that comes directly from the Eynard-Mehta 
theorem. Let 

+£l,2s+£2 (*^! 2r+£i<2s+e 2 (8) 

where -0 is defined in 
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Let w = (w 1: w 2 ) G W, b = (b 1: b 2 ) G B, u x — w 2 + 1, u 2 = (w 2 -w 1 + 2m + 
l)/2, = 6 2 + 1, v 2 = (b 2 -b 1 + 2m + l)/2 and denote 

Sgn(w,b) = (_1)(™™+&i-&2+2)/4 

Define 

/i(w,b) = -sgn(w,b)^ t , 1)Ul (u 2 ,«2) 

and 

2m+l 

/ 2 (w,b) = sgn(w,b) ^ Wl ,2n+i(u2,i - - l)(A~ 1 ) iij ^o,« 1 (j - m- 1,« 2 ) 
where 

/ _ \ 2m+l 

A= (V> ,2n+i(*- "i- !>J -"I -1 ) ) • ( 9 ) 

V / i,j=l 

We then get 



- (-l^-^-'^K^iv^v^m,^) = A(w,b) + / 2 (w,b) 
and we want to prove that 

K- 1 (w,b) = A(w,b) + / 2 (w,b). 
To make the computations simpler, we define T a and C with 
K a (b,w) = -(-l)( fel+b2 - 1 )/ 2 T a (b ! w) 

and 

/>, b) + ; 2 (w, b) = -(-ly^-^c^, b). 

and we will write 

/ J (w,b) = -(-l)(^- 1 )/ 2 /,(w,b). 

Therefore, showing K a .(fi + f 2 ) — I is equivalent to showing T a .C = I. 

We will use the notation that b = b 2 ) and y = (yi,y 2 ) are black 
vertices. We have that 

(r a / i )(b,y) = ^;T a (b,w)/ i (w,y) (10) 
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for % G {1,2} where w ~ y means that w is nearest neighbors to b because 
T a (b, w) = if b and w are not nearest neighbors. We can then write 

(T Q c)(b,w)= ( T «fi)&y)- 

»e{i,2} 



The number of terms on the right hand side of equation ( 10 ) is dependent 
on the location of b and so we split the computation for finding T a C(b,y) 
into the different locations of b. These are given by 

(i) the interior, labeled J, 

(ii) the left hand boundary, labeled L, 

(iii) the bottom boundary, labeled 23, 

(iv) the top boundary but not equal to (2n, 2n — 1), labeled T and 

(v) the special point, (2n, 2n — 1). 

The left hand boundary, L, consists of vertices b = (0,62) where 62 £ 
2Z + 1 and 1 < &2 < 2n — 1. For b G L, we have that b has two neighboring 
white vertices given by b + e\ and b — e^. 

The bottom boundary, 23, consists of vertices b = (61, — 1) where b\ G 2Z 
and 4m + 2 < h\ < 4m + 2n. For b G 23, we have that b has two neighboring 
white vertices given by b + e\ and b + e 2 . 

The top boundary, T, consists of vertices b = (b\,2n— 1) where b\ G 2Z 
and 2n + 2 < b\ < 4m + 2n. For b G T, we have that b has two neighboring 
white vertices given by b — e\ and b — e%. 

For the special point, b = (2n,2n — 1), we have that b has three neigh- 
boring white vertices given by b + b — e2 and b — e\. 

The interior, J, is given by the remaining vertices. For b G J, we have 
that b has four neighboring white vertices given by b ± e r for r G {1,2}. 



In each of the above cases, we evaluate (10). Due to the formulas for fi 
and /2 being rather complicated, we used computer algebra to help with the 
computations. We give the calculation for the first case with full details and 
for the remaining cases, we provide an overview of the main steps. We now 
proceed with checking the above cases. 
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The interior 

Using (|6| and the definition of T a (b, w), we have that for b G J 



T a (b, w) 



±1 if w = b ± e\ 
=pai if w = b ± e 2 
otherwise. 



When b e J and y G Bad, we have that ([lOJ) reads 

T a fi(b, y) = /<(b + ei, y) - /»(b - e u y) + (/,(b - e 2 , y) - / f (b + e 2 , y))ai (11) 
for i e {1,2}. 



We first simplify (11) for i — 1. We can expand out the definition of fi 



in terms of ip. This means that we can rewrite (11) for i = 1 in terms of ip. 
We obtain 



T a / 1 (b,y) = -(-l) 



(bi-b 2 +yi-y2)/4:- l yi+y2 



a ^2/2 + l,&2 
+ ^2/2+1,62 



2m + 1 — yi + y 2 2m — 1 — b x + 6 2 

2 ' 2 

2m + 1 - yi + y 2 2m + 1 - b\ + b 2 
2 ' 2 

2m + 1 - yi + y 2 2m + 1 - 6i + 6 2 
2 ' 2 

2m + 1 - yi + y 2 2m + 3-b 1 + b 2 



(12) 



where ^ is defined in ([8]). We shall evaluate T a /i(b, w) in three cases: b 2 
y 2 ,b 2 > y 2 and b 2 < y 2 . 



For y 2 = b 2 , we only need to consider the last two terms of (12) because 



the first two terms involving ip are zero by <\8h. Using rt8^, we can rewrite (12) 



in terms of ip and hence write each expression as an integral. We find for 
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r„/i(b,y) 



b l+Vl 
I 4 1 



3/2+1,2/2+2 



2m + 1 - yj. + 2/2 2m+l-bi + y 2 



— aip, 



2/2 + 1, J/2+2 



2vri 



2 ' 2 

2m + 1 — 2/i + 2/2 2m + 3 - b x + y 2 % 



0,a 



2vri 



j2i/l 





w 2 

'To, a 

1 if 6i = yi 
otherwise 



W l(bi-Vi) 

a — w 

iiX-i/l-2 



aw 2 



a — w 



dw 



because 2/1 G 2Z. 

As both b and y are black vertices we have that b 2 and y 2 are both odd 
integers. Therefore, the condition that b 2 > 2/2 is equivalent to 62 > 2/2 + 1- 
For 62 > 2/2, all four terms of (12) involving ?// are nonzero and each term can 
be rewritten using if). We find that for b 2 > y 2 



T B / 1 (b,y) = -(-l) 



(61 - 6 2 +2/1 -2/2 ) /4 • yi +V2 



°V'2/2+l,62 

+ ^2/2+1,62 
— ^2/2+1,62+2 



2m + 1 — 2/1 + 2/2 2m — 1 — 61 + 6 2 



2m + 1 - 2/1 + 2/2 2m + 1 - 61 + 6 2 



2 

2m + 1 



(13) 



2/1 + 2/2 2m + 1 - bi + 6 2 



+ aip. 



2/2+1,62+2 



2 2 
2m + 1 - 2/1 + 2/2 2m + 3 - 61 + 6 2 



To evaluate (13), we need to evaluate an expression of the form 

aif)2r,2s+l{xi, X2) + ll>2r,2s+l(Xl, %2 + 1) 
- Ip2r,2s+3{X1,X 2 + 1) + aip 2r!2s+3 (x 1 ,X 2 + 2) 



(14) 



for r, s, Xi, £2 £ Z. We can expand (14) in terms of its integral decomposition 
and combine all the terms under one integral. We obtain 



To, a 



dz 
2~7ri' 



[I + az) 



'\ — a\s-r+l 



1 z 1 (1 + az) az~ 
^ L _| 

z 1-s 



[I + az) 
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In the above equation, the term inside the parenthesis is zero. This means 
we can write 



atp2r,2s+l(xi, X 2 ) + T/>2r,2s+l(zi, ^2 + 1) 
- i>2r,2s+3{xi, X 2 + 1) + a^2r,2s+3(^l, £ 2 + 2) = 



for r < s. 



(15) 



Using the relation in (15), we have that right hand side of (13) is equal to 
zero for b 2 > y 2 which means that T a /i(b, w) = for b 2 > y 2 . 

As b 2 and y 2 are both odd integers, the condition that b 2 < y 2 is equivalent 
to b 2 + 1 < y 2 . We can expand ( 12 ) using the definition of ip given in ^ and 
we find that all the terms of ( 12 ) are equal to zero. Therefore, we have that 
r a /i(b,w) = for b 2 <y 2 . 



We have shown that for b E "3 
T a fi(b,j)-- 



1 ifb = y, 
otherwise. 



For the term T a f 2 (b, y), using (11) we can write T a / 2 (b,y) under one sum. 
We obtain 



T a / 2 (b,y) 



H-^+yi-y 2 , yi+y2 



2m+l 



1> 



l+i/2,2n+l 



2m + 1 - yi + y 2 



'13 



^o,b 2 +i [j-m-1, 



i — m — 1 



2m + 2 - bi + 6 2 



where 



2m + 2 - b x + 6 2 



2m - 1 - 6i + 6 2 



— ^0,62+2 
+ a^o,6 2 +2 



j — m — 1 



2m + 1 - 61 



2m + 1 - 61 + 6 2 



j-m-1, 



2m + 3 - 61 + 6 2 
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As 62 > 0, we can use the relation given in (15) to find that 



A 7 ( ■ 1 2™ + 2-b 1 + b 2 \ 
AjVo,6 2 +i U-m-1, J = 0. 

Therefore, we have 

T a / 2 (b,y) = bGJ,yeB AD . 
To summarize, we have 



T a C(b, Y ) 



The left hand boundary 



1 ifb = y 

otherwise. 



Next, we check T a .C for b on the left hand boundary. For b G £ we have 
that 

1 if w = b + ei 
T a (b, w) = { ai if w = b — e 2 
otherwise. 



For bs£ and y G Bad, using the above equation we find that (10) is given 
by 

Tafifr, y) = /i(b + e 1; y) + f t (b - e 2 , y)oi (16) 

Similar to the interior, we can expand T a /i(b, w) in terms of ip and rewrite 
T a /i(b, w) in terms of an integral using the definition of ip given in ([5]). By a 
computation, we find that for 5 2 = V2, we obtain 

I r e 3niyi/4: w -y 1 /2 

r o /i(b,y) = 7^ dw 



27ri Jr , a w ~ a (17) 
1 if 1/1 = 
otherwise. 

because y\ e 2N, as the integrand has no pole at infinity. For b 2 > y 2 , we 
also find by computation that 

Tahib, y) = — (1 + a 2 ) 

27T1 
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As b 2 > 2/2 (i- e - b 2 > y 2 + 2) and yi > 0, the integrand has no pole at infinity 
or — 1/a, hence 

T o / 1 (b,y) = 0. (18) 

We find that for b 2 < y 2 

r a A(b,y) = (19) 

by using the same reasoning as the case for b in the interior, that is, each 
term in the expansion of T a f\(b, y) in terms of ip is equal to zero. 



Using equations (17), (18) and (19) we have 



T a /i(b,y) 



1 ifb = y 

otherwise 



(20) 



for b G £. Similar to the interior case, using the expansion of T a f 2 (b, y) given 
in (16), we can expand using the definition of / 2 (b,w) to obtain 



T a / 2 (b,y) = (_l)(yi- fe 2-^)/ 4 iyi+^ ^ (A' 1 
2m + 1 - yi + y 2 



,2n+l 

where 

A £ ^o,6 2 +l 



,i-m-l A £ -0o,62+i [j-m-1 



2m + 2 + b 2 



(21) 



j — m — 1 



2m + 2 + b 2 



aipo,b 2 [j-m-1 



2m-l + b 2 



ipo,b 2 +2 (j-m-1 



2m + l + b 2 



We can rewrite the right hand side of the above equation in terms of its 
contour integral using ^ which gives 

1 



2ni 



(1 + a 2 )w- 1+ >- 2m (w - a)(- 3 - fe )/ 2 (l + aw^-Wdw. 



ro,a 



Since 1 < j < m, b 2 > 1 the integrand has no pole at —1/a or infinity, hence 
the above quantity is zero and so the right hand side of (21 ) is equal to zero. 
Combining with (20) gives 

T a C(b, y) = | other J se 

for b G £ and y G Bad- 
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The Bottom boundary 

We now consider the bottom boundary. We have that for b G 23 



T a (b,w) 



1 if w = b + ei 
— ai if w = b + e 2 
otherwise. 



Using the above equation, (10) can be rewritten for b e S and is given by 

T a fi(h, y) = fi(b + ei, y) - ai£(b + e 2 , y). (22) 

We first consider T a /i(b, w) for b 2 — y% = —1. Similar to the analogous 
computation for b in the interior, we can expand the right hand side of ( 22 ) 



in terms of ip and rewrite the expression as an integral. By a computation, 
we find that 



r a /i(b,y) 



2?ri 

1 ifb = y 

otherwise. 



w (!>i-Vi-2)/2 dw 



r ,, 



(23) 



For 62 = — 1 < U2i using the same reasoning as given for b in the interior, we 
have 

T a /i(b,y) = 0. (24) 



Combining equations (|23|) and (|24|), we have 

1 if bi = yi 



r a /i(b,y) 



otherwise 



(25) 



for b G S and y G Bad- For T a f 2 (b, y), using (22) and following the same 



steps for the analogous computation for b in the interior, we have 

2m+l 

1 



T a / 2 (b,y) = (-i)^- b 2+yi~v2)/^+y2 J2 (A' 1 
'2m + l-|/i + y2 



j/2+l,2n+l 



0,2+62 



j - m - 1, 



2m 



z — m — 1 



61 
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where 

2m + 2 - hi + b 2 
2 



A^o,2+6 2 [j-m-1, 



. . 2m + l-6i + 6 2 \ 7 / 2m + 3-6i 
- "00,62+2 ( j — m — 1, + af 0M+2 [ j — m — 1, 



We can rewrite the above equation using the integral definition of if). A 
computation gives 

A,*** (j - rn - 1, 2n, + 2 - 6 ' +62 ) = ± = o 

because b\ > 4m + 2 by the co-ordinates of the bottom boundary. We have 
obtained 

T a / 2 (b,y) = (26) 



1 ifb = y 

otherwise 



for b G 23 and y G Bad- Using (J25J) and (|26|) gives 

T a C(b,y) = 
for b G S and y G Bad- 
The Top Boundary 

We can now consider b = (bi,b 2 ) on the top boundary which means that 
b 2 = In — 1 and bi > In. We have that for be T 

— 1 if w = b — e\ 
T a (b, w) = l ai if w = b — e 2 
otherwise 



We can use the above equation to rewrite (10). We obtain for b G 7 and 
y e Bad 

TJiib, y) = -fi(b - e l5 y) + ai/,(b - e 2 , y) (27) 

We first consider T ai /i(b,y) for 6 2 = y 2 — 2n — 1. By using the analogous 
expansion for b as the interior case in terms of ip an d its integral definition, 
we find that for 6 2 = y 2 = 2n — 1 

T B / 1 (b,y) = 0. (28) 
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For y 2 < b 2 = 2n — 1, we find that 



T a /i(b,y) 



27Ti 

(1 _ aN(l-2n+y 2 )/2 ( )( 2n-l- W )/2 (29) 

wJ dw. 



ro,a 



u ,(2n+l-6i+yi-j, 2 )/2 



For T a / 2 (b,y), using (27) we can follow the analogous computation in the 
interior case, we find that 

2m+l 

T / 2 (b,y) = (-i)(h-b2+yi-V2)/^vi+vu J- . 

7 {2m + l- yi + y 2 . \ - / 2m-& 1 + 6 2 

V ) 2 / 2 +l,2n+l ( g ' 1 ~ m ~ 1 J A ^°- fe 2 I J — — 1, ^ 

(30) 

where 

A 7 /. 2m-b 1 + b 2 \ 7 /. 2m-l-bx + b 2 
At^o,6 2 I J - m - 1, J = a^ ,b 2 ( J ~ m - 1, 

7 /. 2m+l-& 1 + 6 2 
+ -0o,6 2 I J — m — 1, 



2 

1 /• ^-a+i-am-n+ti^^ + aw) 



-dw 

which is found by expanding out the integrands using the integral definition 
of ip given in (J5]) where we have used the fact that b 2 = 2n — 1. Notice that 
we can rewrite the right hand side of the above equation as 

1 f w -2+j-2m-n+6i/2^ + aw )n 

-dw = 



2Wr (l-^) n 

( 2m-b 1 + 2n\ ~ ( 2m + 2 - b x + 2n 

V0,2n+l \J - m - !) n _ a^0,2n+l j - TO - 1, 



2 

(31) 

which follows by using expanding the right hand side of the above equation 
using the integral definition of ip. By definition of the matrix A given in (j9]), 
we have 

m 

^2(A~\ j ip 0j2n+1 (j ~ m - 1, fc) = 5 ijk +m+i- (32) 
i=i 
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Using (32) and (31) in (30) gives 



2m + I — 2/1 + 2/2 2m — b\ + 2n 



y 2 +l,2n+l 



— a ip. 



y 2 +l,2n+l 



2 2 
2m + 1 - yi + y 2 2m + 2 - h + 2n 



2 2 

(_l)(l+6l+2/i-2/2)/4 iy i+y 2 -n /• ft _ o j(l-2n+Ki)/2 ^ + ^ ( 2 „-l-y 2 )/2 



2vri 



(2n+l-bi+yi-y 2 )/2 



Therefore, for y 2 < b 2 using (29) and (33) gives 



r a .(/i + / 2 )(b,y)=0. 
For j/2 = &2 = 2n — 1, using (p3j) we have 



r„/ 2 (b,y) 



2/1 



4 



(^-^-2)72^ 



1 if b x = yi 
otherwise. 



(33) 



(34) 



(35) 



Finally, using (j28j), Q34J) and Q35J) , we have 

1 if 6i = yi 



T a C(b,y) 



otherwise 



for b G T and y G Bad- 



The Special Point 

Finally, we need to consider the special point, i.e. when we choose b = 
(2n, 2n — 1). The special point has three neighboring white vertices and we 

-1 if w= (2n- l,2n-2) 

T( /.,,, o,, -n w) J ~ a[ if w=(2n-l,2n) 
'"< - ' '-^ S ai ifv=(2n + l,2n-2) 

otherwise 



For b = (2n, 2n — 1) and y G Bad, (10) becomes 

Tafiib, y) = -/ f (b — ei, y) + ai/i(b - e 2 , y) - ai/ f (b + e 2 , y). (36) 
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For y 2 = 2n — 1, using the analogous steps for b in the interior, we can write 
using (36) and the integral definition of ip 



(_lW4jn+i/i r aw (2n-l- yi )/2 



27ri 



dw = 0, 



a — w 



(37) 



since yi > 2n + 2, when y 2 = 2n — 1. For y 2 <2n — 1, we can write T a /i(b, w) 
using the integral definition of ip given in ([5]). We obtain 



r„/i(b,y) 



+V2 



(_^(l+J/i-J/2)/4jyi 

2?ri 

(1 - ^)(- 1 - 2 «+^)/ 2 (1 + aw )(2n-l-y 2 )/2 



(38) 



ro.a 



w 



(j/i+l-Jte)/2 



We can now expand out T a / 2 (b, y) using (36 ) and the analogous computations 
given in b in the interior. We find that 



2m+l 



T a / 2 (b,y) = (-i^-fc+w-KO/^i+w ^ (A- 1 ),,,- 
2m + 1 - yi + y 2 



■0 



1+2,2, 2n+l 



where 



z - m - 1 A s ^ ,fe 2 j-m-1 



2m — 61 + 65 



(39) 



AsV'o.fe [j-m-1 



2m — bi + 62 



«^o,fe 2 (j-m-1 



2m - 1 - 61 + 6 2 



+ ^0,62 



2m + 1 - 61 + 62 



2m + 3 - 61 + 6 2 



By setting b = (2n, 2n — 1), we have by a computation that 



A s ^o,2n j-m-1, 



2m - 2n + 2n - 1 



1 

2vri 



w 



-2+j"-2mQ _|_ 



aw) 



^o,2n+i (j - m - l,m + 1) 



-dw 
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where the last equation follows from (|5T) and ((81). We can substitute the above 



equation into (39) and use equation Q32J) in a similar way to the analogous 
computation found in the previous subsection. We find that 

'2m + l-yi + y 2 



T a / 2 (b,y) = (-l) (6l - fe2+ ^ 2) ^ 1+ ^i +J/2) 2n + i 



m + 1 



(_l)(l+M-l«)/4jtfl+i/2 /• w bn-Vl-l)/2(l _|_ att ,)(2"-l-!/2)/2 



27Ti 



r ,, 



(i-S) 



a x(2n+l-y 2 )/2 



(40) 



We have that for y 2 = 2n — 1, (|40|) is equal to 
T a / 2 (b,y) 



1 

2txi 



(— 1) 2/1 / 4 i™+^ 1 ^^-i+i/i / 2 ^. 



w 



ro,a 



1 if 2/1 = 2n 
otherwise. 



for b = (2n,2n - 1). 

Therefore, we have for b = (2n, 2n — 1) 

T a .C(b,y) = ^ T a /i(b,y) 

i6{l,2} 



1 ifb = y 

otherwise 



by using equations (37), (38) and (40) for b = (2n, 2n — 1) and y G Bad- 



7 Proof of the formula for the L-kernel 

In this section, we derive the L kernel from the inverse Kasteleyn matrix. 
Let bj = (£j,Tjj), 1 < j < I, be the positions of black vertices in Kasteleyn 
coordinates. We want to show that 

P [L-particles at &i, . . . , bi] = det (h n ^ p 77^ £j, Vj)\<ij t <i (41) 

Note that we have an L-particle at a black vertex of b if and only if a 
dimer (domino) covers the edges (b,b + f\) or (b,b + / 2 ) where fi = t\ and 
j'2 = — e 2 . Hence, by Theorem 272] and the linearity of the determinant in its 
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rows, we have that 

2 i 

P [L-particles at 61, ... , 6,] = ^ fj #„(&(, h + f n ) det (K" 1 ^ + frM)^.^ 

ri,—.n J=l 

= det ( ^ iT a (6,, 6 4 + /JiC 1 ^ + /r t , ^ 

V^=i / i<i,j<i 

Recall that, from ([6]), we have that 

K a (b l M + h) = (-ir^ +Th+1)/2 

and 

K a {hX + h) = (-l)-^ +l)/2 ai. 
If we use the formula ^ for the inverse Kasteleyn matrix, we see that 

2 

(-l)(*-6)/*-(m-&)/4 J- ifa ( 6lf ^ + / r )K~ 1 (6i + / r ,6 2 ) 

r=l 

= Kn,p (7/2 + 1, 2 ' 771 + ' 2 

- aK n:P [772 + I, 5771, 

where M = 2m + 1 and we have used (— = — 1 since 771 is odd. Thus, 



to prove (41) it suffices to show that 



L„, p (4i, 771; 772) = JK„ iP (772 + I, ; r/i + 2, 



- aK„ iP (772 + I, ; 771, 2 

Now, with £j G {0, 1} we have from ^ we have 

+ (-l) x -^(2r + e 1 ,x;2 S + e 2 ,7/) 

- (-1)-* (((I - K) 2r l n+1 a^ 2s+£2 )(k), b^, 2r+£l (k))f (2m+l) 
2 

■.= J2 R nU 2r + £ ^^ 2s + £ ^y)- 

8=0 
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Define 

T^ m {£l,Vi;&,V2) = Rn!m [V2 + h ^ ^1+2, 

- a K; m \ V2 + i, g ; 7/1 ' 2 J ' 

(42) 

We have to show that 

2 

^ Ti%(2r + £l , x; 2s + e 2 , y) = Ln,„(£i> W 6, (43) 

i=0 

Since r\\ and 772 are odd we can write 771 = 2s — 1, 772 = 2r — 1 so that 
772 + 1 = 2r and 7/1 = 2(s - 1) + 1 and 7/1 + 2 = 2s + 1. Now by (J5j, 



i 2s+K2r 



( 1 + a2 )' r r fe-^l)/2+(6-6)/2^ 



-2 



2ttz 7 ro a (1 - a/^-'+i z (44) 

al28-K2r /" (1 + a^) S ~ 1 ~ r ^ fa2 _ ;?l) /2 + fa- £2 )/2 + l^ 

2tu y ro a (1 - a/zy-r z ' 

If 771 > 772, i.e. s > r + 1, then 2s + 1 > 2s — 1 > 2r and the expression in 
the right hand side of (44) is equal to zero. If 771 = 772, i.e. r = s, we get 



a 



-dz = — _ . / — dz 



2ni J To 1 + az 2ixi J Vq 1 + az 

because of the change of orientation from moving the contour Tq^ to the 
contour Y_y a . Thus, 

rm / s l + a 2 f (1 + az)*™-**)/*- 1 

JI Cl<€2 - U -T ? i=r 72 k a ) 
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If £1 > £2 and 771 < 772 the integral in (45) is equal to since the integrand is 
analytic inside r , a - Thus, the first term in (45) equals 

1 + a 2 



+ (6 



"(1 l»?i>r/2 ~^-vi=m 



6)/2 d,2 

r0ja (* - a)(*-*0/»+i ~ 

27rz y roa ( z _ )(m-^)/2+i 2 



<6 



The last integral is equal to zero if rji > 772 and hence the last expression 
equals 



1 + a 2 



2tu 



2 («i-6)/2 



-1/a 



(1 + az)(z — a) 



dz 



L «i<6 



1 + a 2 
2tu 



r 



(l + az)fa-^-i (£1 _ £2)/ 
(2 _ a ) (771-^/2+1 z 



The first term in this expression equals 



-^1=772 -"-£i<6 



Ki-6)/2 



so combining this result with (45) we find 

T n°i(6,^i;6,^2) 



-1 



1 + a 2 
2vri 



(1 + a^)^- 7 ' 2 )/ 2 - 1 Ki . 



-6)/2 



(46) 



, r0ia (z-o)(^-*)/^ 
Next, using ([5]) and inserting it into (42) we get 
7^,1(6, W 6, 772) 

dv t; n-6/2 (! + an )(m-l)/2( n _ a )n-(^+l)/2(47) 



1 + a 2 

(27ri) 2 



dv 



, o o u - u u 11 -^/ 2 (1 + at>) (7?2+i:)/2 (i> - a)™-^ 2 - 1 )/ 2 
after a short computation. Finally, by (42) and ^ we see that 

= ((I - K)2^ +1 (a_( m _£ 1+ j W -)/2 ! 2 s +l + aa_( m _£ 1+ M-2)/2,2s-l)(fc), 



"(»72 



-b+M)/2M K )/ >2 m+l { l > 
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Now, using ([5]) a computation gives 

(_-Q-(m-«i+M-2)/2 ( a _ (r?i _ 6+M)/22s+1 (/ c ) + aa_ (7?1 _ a+M -2)/2 > 2 S -i(A;)) 

l) fe (l + a 2 ) f , f dv V (ti-Vi-l)/2(l +av yVi-l)/2(l_ a / v \n-( m + l)/2 



(2tu) 2 7 roa Jr , a , u u-v u k + l (1 + ov)»(« - a)«+ x 

;i + a 2 )% fll (A;). 

(48) 



Note that by moving the t> -contour in (48) inside the w-contour we get the 
expression in (|4]). Similarly, by 

&-fe-6 + M)/2,2,(fc)(-l) fe - ?2+M)/2 

^-l) fe [ d [ dv v k (l + av) n (w -a) n+1 



(2m) 2 J Toa Jr 0au v ~ uu^-^ +1 )l 2 (1 + av)to>+ 1 V 2 (l - a/v) n -^- 1 )/ 2 

Again by moving the f-contour inside the w-contour we get the expression 
in Q. Thus, 

TSU^Vi,^) = ~ <((I " K)^ +1 A ium )(k),B^ V2 (k)) >2m+1 (1 + a 2 )- 

(49) 



If we use (46), (47) and (49) we obtain (43) which is what we wanted to 
prove. 
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